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Abstract 

Using a renormalization group method, we calculate 800 high-temperature 
coefficients of the magnetic susceptibility of the hierarchical Ising model. The 
conventional quantities obtained from differences of ratios of coefficients show 
unexpected smooth oscillations with a period growing logarithmically and can 
be fitted assuming corrections to the scaling laws with complex exponents. 
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The renormalization group method has enhanced considerably our understanding of 
elementary processes and critical phenomena. In particular, it has allowed the computation 
of the critical exponents of lattice models in various dimensions. On the other hand, the 
critical exponents can be estimated from the analysis of high-temperature series 0. Showing 
that the two methods give precisely the same answers is a challenging problem ||. More 
generally, much could be gained if we could combine these two approaches, in particular, in 
the context of lattice gauge theories. 

As far as the numerical values of the critical exponents are concerned, there are two 
difficulties. The first one |4j] is that one needs much longer high-temperature series than the 
ones available || (which do not go beyond order 25 in most of the cases) in order to make 
precise estimates. The second is that the practical implementation of the renormalization 
group usually requires projections into a manageable subset of parameters characterizing 
the interactions. It is nevertheless possible to design lattice models ||, called hierarchical 
models, which can be seen as approximate version of nearest neighbor models, for which such 
projections are unnecessary. For the hierarchical models, the renormalization group transfor- 
mation reduces to a recursion formula which is a simple integral equation involving only the 
local measure. This simplicity allows one to control rigorously |7[ the renormalization group 
transformation and to obtain accurate estimates of the critical exponents |§ . As explained 
in the next paragraph, the recursion formula also allows to calculate the high-temperature 
expansion to very high order. Consequently, the hierarchical model is well suited to study 
the questions addressed above. 

In this letter, we report the results of a large scale calculation which performs the high- 
temperature expansion of the magnetic susceptibility of the hierarchical Ising model up to 
order 800. The method of calculation uses directly the renormalization group transformation 
with a rescaling of the spin variable appropriate to the study of the high-temperature fixed 
point. This method has been presented in Ref. |9| and checked using results obtained with 
conventional graphical methods [ItJ. For the sake of briefness, we shall follow exactly the 
set-up and the notations of Refs. PHTTI where the basic facts concerning the hierarchical 



model were reviewed. The hierarchical models considered here have 2 n sites, and the free 
parameter c which controls the strength of the interactions is set equal to 2 l ~ 2 l D in order 
to approximate a nearest neighbor model in D-dimensions. We define the finite volume 
magnetic susceptibility as 

Xn((3) = 1 + h, n /3 + b 2 ^ 2 + ... (1) 

In most of the calculations presented below, we have used n=100, which corresponds to a 
number of sites larger than 10 30 . From a mathematical point of view, the calculation of Xn 
amounts to repeating n times the Fourier transform of the recursion formula |J 

Rn + i(k) = C n+l eM-\^l) n+1 ^){Rn{^)f , (2) 

expanded to the desired order in ft, with the initial condition R Q = cos(fc). The constant C n+ i 
can be adjusted in such a way that the Taylor expansion of R n+ i(k) reads 1 — (l/2)k 2 x n +i + ---- 
This calculation has been implemented with a C program and ran for 6 weeks on a DEC- 
alpha 3000 in order to obtain 800 coefficients for D = 3. 

For the discussion which follows, it is crucial to estimate precisely the errors made in 
the calculation of the coefficients. There are two sources of errors: the numerical round-offs 
and the finite number of sites. We claim that with 2 100 sites and 3 < D < 4, the finite 
volume effects are several order of magnitude smaller than the round-off errors. From Eq. 
(2), one sees that the leading volume dependence will decay like (c/2) n . This observation 
can be substantiated by using exact results at finite volume |H| for low order coefficients, 



or by displaying the values of higher order coefficients at successive iterations as in Figure 
1 of Ref. ||. In both cases, we observed that the (c/2) n law worked remarkably well. For 
the main calculation presented below, we have used c = 2i (i.e. D = 3) and n = 100, 
which gives volume effects of the order of 10~ 20 . On the other hand, the round-off errors are 
expected to grow like the square root of the number of arithmetical operations. In Ref. [|9] , 
we have estimated that this number was approximately n.m 2 for a calculation up to order m 
in the high-temperature expansion with 2" sites. Putting this together, we estimated that 



for n = 100, the error on the m-th coefficient will be of order m x 10 -16 . We have verified this 
approximate law by calculating the coefficients using a rescaled temperature and undoing 
this rescaling after the calculation. We have chosen the rescaling factor to be 0.8482 and 
the rescaled critical temperature is then approximately 1. This prevents the appearance 
of small numbers in the calculation. If all the calculations could be performed exactly, 
we would obtain the same results as with the original method. However, for calculations 
with finite precision, the two calculations have independent round-off errors. Comparing 
the results obtained with the two methods for the coefficients up to order 200 shows that 
the numerical fluctuations of b m grow approximately like m x 10~ 16 . More conservatively, 
we can say that the numerical errors are bounded by m x 10~ 15 . We conclude that for the 
calculations reported below, the errors on the coefficients are dominated by the numerical 
round-offs and we estimate that they do not exceed 10~ 12 . 

In order to estimate 7, we used standard methods described in Refs. [@,|J. For the sake 
of definiteness, we recall a few definitions. First, we define r m = b m /b m _i, the ratio of two 
successive coefficients. We then define the normalized slope S m and the extrapolated slope 
S m as 

S m = -m(m - l)(r m - r m _i)/(mr m - (m - l)r m _i) ; 

(3) 

S m = mS m - (m - l)S m -! . 

The extrapolated slope, which is free of order n~ l corrections ||, is displayed in Fig.l for 
m < 200. For comparison, we have also displayed the results for D = 3.5 and 4. A surprising 
feature is the clear appearance of large oscillations for D = 3. When D is increased, the 
amplitude of these oscillations diminishes. They are still present at D = 4 and can be seen 
better by plotting S m+ \ — S m . One important point of this letter is to establish that these 
oscillations are not due to the errors discussed above. As a consequence of the multiplications 
by m appearing in the definition of the the extrapolated slope, the errors are amplified by a 
factor which can be as large as 10 5 for m near 100 and 10 7 for m near 500. However, even 
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when multiplied by such a large factors our most conservative estimate of the numerical 
errors gives errors on the extrapolated slope which are several order of magnitude smaller 
than the amplitude of the oscillations. We have made independent checks of this statement 
for D = 3 by calculating directly the extrapolated slope for n =100 and 200 and by using an 
intermediate temperature rescaling as explained above. The smoothness of the oscillations 
appears clearly in Fig. 2, where S m is displayed for 50 < m < 800. This smoothness rules 
out large numerical fluctuations. In conclusion, we have established that the oscillations in 
the extrapolated slope are a genuine feature of the model considered. Fig. 2 also shows 
that the extrema are not equally spaced. Instead, the location of one extremum can be 
approximately found by multiplying the location of the previous extremum by 1.19. In 
other words, the extrema of Fig. 2 would look equally spaced if the abscissa variable had 
been ln(m) instead of m. This of course suggests the use of a complex exponents since 
Re(m la ) = cos(aln(m)). 

In the conventional description [^] of the renormalization group flow near a fixed point 
with only one eigenvalue Ai > 1, one expects that the magnetic susceptibility can be ex- 
pressed as 



with A = |Zn(A2)|//^(Ai) and A2 being the largest of the remaining eigenvalues. It is usually 



If A is real, there is no room for the oscillations in this description. Nevertheless, the 
fact that the period of oscillation increases logarithmically with m suggests that one could 
modify slightly Eq. (5) by allowing B and A to be complex and selecting the real part of 
the modified expression. This introduces two new parameters and we have chosen to use 
the following modified parametrization of the extrapolated slope: 



x =(f3 c -(3)-y(A + A 1 ((3 c -(3) A + ....) , 



(4) 



assumed that these eigenvalues are real. This implies Hj that 




(5) 
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This parametric expression allows good quality fits for m large enough. For instance, a 
least square fit for the m > 300 data, yields p = 1.412, m = 512, p = 0.67, 7 = 1.310 
and K = 2.53. The fit is displayed on Fig. 2. More accurate results could presumably 
be obtained if we had a consistent description of the oscillations involving definite relations 
among the parameters of Eq. (6). The value of 7 is in good agreement with the result [7|,|3| 
obtained with the e-expansion, namely 1.300. The value of p is not far from |Zn(A2)|//^(Ai) 
which is approximately 0.46 according to Ref. ff[]||. 

We have considered two possible explanations of the oscillatory behavior. The first one 
is that one could replace A2 by a couple of complex conjugated eigenvalues. This possibility 
is not realized in any perturbative calculation we know. For instance, the gaussian spectrum 
is real and its largest eigenvalues are widely separated. One could imagine that when D 
is continuously evolved from 4 to 3 for the hierarchical model, two real eigenvalues merge 
into each other and subsequently evolve as complex conjugate of each other. However, 
the fact that the oscillations persists at D — 4, as indicated in Fig. 1, goes against such 
an explanation. Conformal theories in two dimensions provide examples of calculations of 
the renormalization group eigenvalues as the eigenvalues of the matrix of derivative of the 
beta functions. For low order calculations, this matrix is symmetric |T3[ which implies real 



eigenvalues. A more attractive possibility, is that the susceptibility of the hierarchical model 
satisfies a renormalization group equation of the type discussed in section II of Ref. |14| . This 

i2TV 

equation allows corrections to the scaling law of the type (l + (Ai(/3 c — (3) ln( ~ xs > +c.c) + ..), in the 
simplified case where only one eigenvalue A is considered. Considering the values obtained 
with the fit we see that p = 1.412 is close to the the value of the largest eigenvalue 
Ai=1.427. This result suggests that one should try to find an equation for the susceptibility 



of the hierarchical model related to the one discussed in Ref. [14]. 

In conclusion, we have shown that a calculational method of the high-temperature ex- 
pansion based on the renormalization group method can be a very powerful tool when the 
hierarchical approximation is used. It would be worth trying to improve this method beyond 
this approximation. Our analysis of the magnetic susceptibility has shown that unexpected 
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oscillations appear in the extrapolated slope. A detailed understanding of these oscillations 
is required in order to allow a precise comparison between the results obtained from the 
high-temperature expansion and the e-expansion . 

One of us (Y.M.), stayed at the Aspen Center for Physics during the last stage of this 
work and benefited from stimulating conversations with the participants, especially with N. 
Warner. 
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Figure Captions 

Fig. 1: The extrapolated slope S m for m < 200 and D — 3, 3.5 and 4. 

Fig. 2: The dots are the extrapolated slope S m for 50 < m < 800 and D = 3. The continuous 

curve is the fit described in the text. 
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